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We have analyzed the theory of the baryon magnetic moments in the (approximate) 
QCD setting suggested by Brambilla et al. By modifying their derivation of the qqq 
interaction and wave equation for baryons, we derived expressions for the baryon 
moments in terms of the underlying quark moments, including the first corrections 
associated with the binding of the quarks in baryons. The results, which hold in 
the "quenched approximation" in which the contributions of virtual quark pairs are 
neglected, fail to describe the measured moments, with typical errors on the order of 
7%. We conclude from our analysis that the quenched approximation is at fault, and 
that the baryon moments give a sensitive test of that standard approximation in the 
lattice and analytic approaches to QCD. 

1 Introduction 

The magnetic moments of the stable baryons are very well measured, and their 
theoretical calculation gives a sensitive test of our understanding of baryon 
structure in quantum chromodynamics (QCD). Although the basic pattern and 
approximate magnitudes of the moments can be explained using the nonrela- 
tivistic quark model, the deviation of the moments from the quark-model pat- 
tern has not been explained dynamically despite many attempts. We discuss 
that problem here from the point of view of nonperturbative QCD. In particular, 
we have analyzed the theory of the moments in the approximate QCD setting 
developed in the work of Brambilla et al. jl| , who used a Wilson-loop approach 
to study the qqq bound state problem. By modifying their derivation of the qqq 
interaction and wave equation, we have derived expressions for the baryon mo- 
ments in terms of the underlying quark moments, including the first corrections 
associated with the binding of the quarks in baryons. Our results hold in the 
"quenched" approximation in which there are no internal quark loops embedded 
in the world sheet swept out by the Wilson lines joining the valence quarks, and 
no pairs associated with the valence lines. The corrections to the moments are 
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of relative order (V)/m q , where V is a typical component of the binding po- 
tential and m q is a constituent quark mass, and are potentially large enough to 
explain the deviations of the measured moments from the quark-model values. 

To test the theory, we have constructed variational wave functions for the 
baryons using the interactions derived by Brambilla et al. jjj , including all spin 

and orbital configurations possible for J p = ^ + and internal and total orbital 
angular momenta L < 2, and used them to calculate the moments. The effects 
of excited orbital contributions to the moments are negligible, as expected. The 
new contributions tend to cancel, and the remnants do not have the correct 
pattern to explain the discrepancies between theory and experiment. 

The approximations underlying the Wilson-line analysis are known. Since of 
these, only the quenched approximation has a direct effect on the spin depen- 
dent terms with which we are dealing, we conclude that the use of the quenched 
approximation is responsible for the deviations of theory from experiment, and 
that the moment problem provides a sensitive test of this standard approxima- 
tion in lattice and analytic QCD. 

In a world-sheet picture, the inclusion of internal quark loops which describe 
meson emission and absorption by the baryon would allow the introduction of 
internal orbital angular momentum and spin, and would affect the moments. 
We are now investigating meson-loop contributions using approximations sug- 
gested by the world-sheet picture and chiral pertubation theory. In the following 
sections, we sketch our derivations and calculations, and justify our conclusions. 
A more detailed discussion will be given elsewhere. 



2 Baryon moments in QCD 
2.1 The problem 

The simple, nonrelativistic quark model gives a qualitatively good description 
of the baryon moments. Under the assumption that each baryon is composed of 
three valence or constituent quarks in a state with all internal orbital angular 
momenta equal to zero, the moments are given by expectation values of the spin 
moment operators 

Mqm = ^2 ^"q' CO 
q 

leading to the standard expressions 

1 e 

A*p = o( 4 ^ ~ Md), ■ ■ ■ , Vq^TT 3 -- (2) 

The masses in the quark moments are clearly effective masses, and can be treated 
as free parameters in attempting to fit the data. A least-squares fit to the 
measured octet moments, taken with equal weights, is shown in Table 1. The 
moments are given in nuclear magnetons (nm). The pattern of the signs of 
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the quark model moments agrees with observation, while the root-mean-square 
deviation of theory from experiment is 0.14 nuclear magnetons, or about 9% 
of the average magnitudes of the moments. Agreement at this level can be 
regarded as an outstanding success of the quark model, but the deviations also 
give a very sensitive test of baryon structure: there is presently no completely 
successful first-principles theory of the moments. 

Table 1: Quark model fit to the magnetic moments in the baryon octet. All 
moments are given in nuclear magnetons. 



Baryon 


Experiment 


Quark Model 




P 


2.793 


2.728 


0.065 


n 


-1.913 


-1.819 


-0.094 


£+ 


2.458 ±0.010 


2.639 


-0.181 


£- 


-1.160 ±0.025 


-0.999 


-0.161 


|E°^A 7 | 


1.61 ±0.08 


1.575 


-0.03 


A 


-0.613 ±0.004 


-0.642 


0.029 


CO 


-1.250 ±0.014 


-1.462 


0.212 




-0.651 ±0.025 


-0.553 


-0.098 



2.2 The Wilson-loop approach to the baryon moments 

Our approach to the baryon moment problem is based on the work of Brambilla 
et al. Q, who derived the interaction potential and wave equation for the 
valence quarks in a baryon from QCD using a Wilson-line construction. The 
basic idea is to construct a Green's function for the propagation of a gauge- 
invariant combination of quarks joined by path ordered Wilson-line factors 

U = Pexpig J A g ■ dx, (3) 

where A g is the color gauge field. The Wilson lines sweep out a three-sheeted 
world sheet of the form shown in Fig. |l| as the quarks move from their initial to 
their final configurations. 

By making an expansion in powers of l/m q and considering only forward 
propagation of the quarks in time using the Foldy-Wouthuysen approximation, 
Brambilla et al. are able to derive a Hamiltonian and Schrodinger equation 
for the quarks, with an interaction which involves an average over the gauge 
field. That average is performed using the minimal surface approximation in 
which fluctuations in the world sheet are ignored, and the geometry is chosen 
to minimize the total area of the world sheet per Wilson, subject to the motion 
of the quarks. The short-distance QCD interactions are taken into account 
explicitly. 
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world sheet 




Wilson line 

Figure 1: World sheet picture for the structure of a baryon. 



The result of this construction is an effective Hamiltonian to be used in a 
semirelativistic Schrodinger equation H^f = e\&, 

H = Y) \fpi + ™Z + cx[r x + r 2 +r 3 ) - V \ ^ 

- O Tij 

% i<j J 



-^- Sl ■ (n x Pl ) + ^ • [(r 12 x Pl )^- + (r 13 x Pl )^ 



2 1 a s 2 1 a s 
3-si • ri2 x p 2 - r si • r i3 x p 3 H . (4) 

jmim2 fj 2 d mim3 rf 3 

Here = Xj — Xj is the separation of quarks i and j, is the distance of quark 
i from point at which the sum r*i + r 2 + r 3 is minimized, and and Sj are the 
momentum and spin operators for quark i. The parameter <r is a "string tension" 
which specifies the strength of the long range confining interaction, and a s is 
the strong coupling. The terms hidden in the ellipsis include tensor and spin- 
spin interactions which will not play a role in the analysis of corrections to the 
moment operator, and the terms that result from permutations of the particle 
labels. The full Hamiltonian is given in H. This Hamiltonian, including the 
terms omitted here, gives a good description of the baryon spectrum as shown 
by Carlson, Kogut, and Pandharipande Q and by Capstick and Isgur who 
proposed it on the basis of reasonable physical arguments, but did not give 
formal derivations from QCD. 

The presence of the quark momenta p^ in the Thomas-type spin-dependent 
interactions in Eq. (^) suggests that new contributions to the magnetic moment 
operator could arise in a complete theory through the minimal substitution 

Pi > Pi — BiA. em {Xi) , (5) 



4 



with A cm (xi) the electromagnetic vector potential associated with an external 
magnetic field. However, this point is obscured in the work of Brambilla at 
al. by the transformations they make to express the equations for the Green's 
function in Wilson-loop form. We have therefore repeated their derivation of 
the valence-quark Hamiltonian in Eq. ( ||), replacing the SU(3) C color gauge field 
gA g by the full SU(3) c xU(l) em gauge field gA g + e q A cm at the beginning. By 
then reorganizing the calculation of the three-quark Hamiltonian to keep e q A em 
explicit throughout, and then expanding to first order in A cm with 



A om = — B x x g , B = constant, (6) 



1 

2 

we can identify the modified magnetic moment operator through the relation 



AH = -n-B. (7) 

The new moment operator, \x = /xq M + Afi, involves the leading corrections 
to the quark-model operator associated with the binding interactions. Afi can, 
in fact, be read off from the terms in Eq. (|]) which depend on both the quark 
spins and momenta by making the minimal substitution in Eq. (j^). For example, 
the term which depends on Si ■ ri 2 x pi gives an extra contribution 

xi x (si x ri 2 )-g- (8) 



6mf 



to fi 1 . There are also possible orbital contributions to the moments because 
the Hamiltonian mixes states with nonzero orbital angular momenta with the 
ground state. These have the standard form to the accuracy we need. 
The baryon moments are now given by expressions of the form 
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i=l i=l 

where we have quantized along B, taken along the z axis. The expectation 
values are to be calculated in the baryon ground states. The correction terms 
Si from the new operators are given for the L — baryons other than the A by 

3e + ei e 3 e 2 A + Ai . 

Si = — + , i = l,2 

Zmi ei rri3 Lrti\ 

5s = iL + ^iiL__^ ; (10 ) 

m 3 e 3 mi m 3 

where the e's and A's are ground state radial matrix elements, 

,2 a s 2a s r 2 3-z 2 , 2 a s r 3 i • z 3 

12 6r 23 6r 3i 
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The identical quarks in these baryons arc labelled 1 and 2, the unlike quark, 
3. In writing these results, we have made the approximation r% + r 2 + r% w 
hi r i2 +^23 +J"3i), known to be reasonably accurate for the ground state baryons 
[pj, and used the corresponding Thomas spin interaction. The result for the A 
is similar. 

2.3 Tests of the model 

Rough estimates of the matrix elements above suggest that til mi w 0.05 and 
that Aj/mj rs 0.3 for the light-quark masses used in Refs. 2 and 3. Binding 
effects are therefore potentially large, and are different for different quarks and 
baryons. To obtain a quantitative test of these effects and reliable estimates 
of the orbital contributions to the moments, expected to be small, we have 
performed a detailed analysis of the baryon wavefunctions using the Hamiltonian 
given in We use Jacobi-type internal coordinates p, A, with p = xi — Xa, 
and A = R12 — X3, where Ry is the coordinate of the center of mass of quarks 
i and j. The most general j = \ + baryon wave function for quarks 1 and 2 
identical and orbital angular momenta L p , L\, L < 2 has the form 

-0±, m = [ a %l) a l + ib Tpb (ci - <y-i) ■ p x A + c o -0 c t\i{p) 

+d ip d t 12 (\) + e ip d Oi - o- 2 ) • 0-3 p ■ A ] xf" =1 . (12) 

2 ' m 

Here t±2 is the usual tensor operator 

t 12 (x) = 'dcri ■ xer 2 • x - cr 1 ■ er 2 x 2 , (13) 

and xf 12=1 is the standard three-particle spinor for Su = 1, j = 5, is = m. The 

scalar functions ^ = t/Ji(p 2 , A 2 ) are normalized together with the accompanying 
spin operators. The constant coefficients ao, ■ ■ ■ , eo, normalized to unity, give 
the fractions of the various component states in With this form of the wave 
function, we can use trace methods for the spins to calculate such quantities as 
(H) and (p). 

We have made variational calculations of the energies and approximate wave 
functions of the ground state baryons and their first excited states using the 
Hamiltonian in [0. The information on the excited states allows us to calcu- 
late the coefficients 60, . . . , eo perturbatively. While the Hamiltonian does mix 
orbitally excited states into the L p = L\ = L = quark-model ground state, 
the coefficients are very small, ranging from essentially zero to about 0.02 de- 
pending on the baryon. Because these coefficients only appear quadratically in 
the baryon moments, the orbital contributions to the moments are completely 
negligible. 

The radial matrix elements €i and Aj defined in Eq. ( |fl] ) are significant, with 
the e's ranging from 10 to 20 MeV, and the A's from 40 to 70 MeV. A new fit 
to the moment data using the expressions in Eqs. (O) and (Eoj) , with the quark 
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masses allowed to vary, gives a small improvement in the fitted moments, with 
the root-mean-square deviation from the measured moments decreasing from 
0.14nm for the quark model to O.lOnm. A secondary effect of the corrections 
is to change the fitted quark masses or moments significantly, with the effective 
quark masses decreasing, or the moments increasing. 

We have concluded from this exercise first, that the finer details of baryon 
structure are not yet described correctly in the present QCD-bascd model, and 
second, that the baryon magnetic moments provide a very sensitive check on 
the theory. 

3 Possible improvements in the theory 

We begin the discussion of possible improvements in the theory of the moments 
by recalling the approximations used in the construction given by Brambilla et 
al. jjj and in our derivation of the moment operators: 

(i) the l/m q expansion; 

(ii) the minimal surface approximation; 

(iii) forward propagation of the quarks in time; 

(iv) and the quenched approximation. 

Of these approximations, only (iii) and (iv) are likely to affect the moments 
significantly. 

The l/m q expansion is to be interpreted as an expansion in constituent 
quark masses. It can be resummed in the kinetic terms in the Hamiltonian as in 
Eq. (||), and leads elsewhere to the appearance of effective inverse masses 1/mj 
which need not be the same as the kinematic masses, but represent averages of 
quantities such as = 1/ y/pf + fnf. Since there is no new spin dependence 
involved, and we have treated the quark masses as free parameters in fitting the 
moments, relativistic corrections of this sort are unlikely to change our results. 

The average over the color gauge field A g in the minimal surface approx- 
imation neglects fluctuations, and minimizes the surface energy of the world 
sheet to obtain the approximate Hamiltonian. Since the color fields do not 
carry charge, they do not contribute internal currents in the baryon, and an 
improved treatment of the averaging would presumably not change the moment 
operator directly. While it could change the functional form of the Hamiltonian 
somewhat, the quantitative success of the model in describing baryon spectra 
suggests that the changes would not be large, and their effect on the moments 
through the e and A parameters would be minimal. 

In contrast, the remaining two approximations affect the moment operator 
directly. Internal quark loops can contribute circulating currents in the baryon. 
These are omitted in the quenched approximation. In addition, the forward 
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propagation of the quarks in time inherent in the use of the Foldy-Wouthuysen 
approximation eliminates quark pair effects connected with the valence lines in 
Fig. 1. We believe that this general quenched picture underlies the difficulties 
with the model, and that it will be necessary to include pair effects to obtain a 
precise dynamical description of the moments. 

One effect of internal quark loops can be seen in Fig. |. In this fi gure, we 
suppose that a quark loop is embedded in the minimal world sheet of the baryon. 
The effect is to give a meson state and a new baryon in a world-sheet picture of 
meson emission and absorption. An external magnetic field interacting with the 



meson 




new baryon 



Figure 2: Quark loops embedded in the world sheet give meson-baryon states. 

system will see a meson current as well as the baryon current, and the moment 
will be modified. Meson currents were invoked in the past in attempts to explain 
the full anomalous moments of the nucleons. Here we are only concerned with 
presumably small corrections to the quark-model moments. 

One approach to the calculation of meson loop effects is through chiral per- 
turbation theory. This has a long history, and has been studied recently in 
the context of baryon moments by a number of authors 0, |H[ ^| . The relevant 
diagrams for the interaction of the baryons with the pseudo Goldstone bosons 
of the chiral theory are shown in Fig. [| The work of Jenkins et al. JiJ uses 




Figure 3: Goldstone-boson diagrams which contribution to the magnetic mo- 
ments. 

heavy baryon chiral perturbation theory J7J, and studies the effect on the mo- 
ments of the "nonanalytic terms" in the symmetry breaking parameter m s . The 
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(a) (b) 



Figure 4: (a) Appearance of Bethe-Salpeter wave functions in the diagrams, (b) 
Schematic reminder of the distributed nature of the meson-baryon vertex. 

analytic terms in m s are ambiguous because the inevitable appearance of new 
couplings at each order in the chiral expansion, and are ignored. The results 
of Jenkins et al. are not encouraging as they stand, but we have found errors 
in some of the coupling factors given in the published paper. Luty et al. || 
concentrate on a simultaneous expansion in 1/N C and m s , and obtain interest- 
ing sum rules for the moments but little dynamical information. Finally, Bos et 
al. H consider the moments from the point of view of flavor SU(3) breaking in 
the baryon octet using a different chiral counting than that usually used, and 
obtain a very successful parametrization for the moments. However, this model 
is again nondynamical, and has seven parameters to describe eight measured 
moments. In fact, a fundamental problem with the chiral expansion from our 
perspective is that it simply parametrizes the moments with an expansion con- 
sistent with QCD, but does not control the higher order terms in what appears 
to be a slowly convergent series. 

We are presently investigating meson loop effects using a somewhat different 
approach suggested by the world-sheet picture. The baryon appears in this 
picture as an extended object which must absorb the recoil momentum in the 
emission of a meson. We would therefore expect wave function effects (form 
factors) to be important for high meson momenta, and to supply a natural 
cutoff for loop graphs. The wave function appears naturally when the process 
is viewed using "old fashioned" instead of Feynman perturbation theory, as 
indicated in Fig. energy denominators and vertex functions combine to give 
the B' M component of the wave function for a baryon B, or, with different 
time ordering, the BM component of the wave function of B' . Since the wave 
functions are expected to be fairly soft, with characteristic momenta below the 
chiral cutoff of ~ 1 GeV, they can be expected to combine a number of terms in 
the chiral expansion in a way that is dynamically accessible through approximate 
models derived from QCD such as that of Brambilla et al. jl| . The non-point 
character of the meson-baryon vertex in spacetime indicated in Fig. ^ contrasts 
sharply with the point vertex used in chiral perturbation theory, Fig. 0, and is 
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also likely to play a role. We have reached the same conclusions by studying the 
exact sideways dispersion relations for the moments given by Bincer |8). The 
problem there is in extracting the quark-model moments. 

We note finally that a different aspect of symmetry breaking, the suppression 
of strange-quark loops through mass effects, appears as a natural possibility in 
a world-sheet picture. Note that this is not the same as the suppression of kaon 
loops considered by other authors; see, e.g., jl and the references given there. 

4 Conclusions 

On the basis of the work sketched above, we have concluded that the magnetic 
moments of the baryons give a sensitive test of baryon structure and of approx- 
imations in QCD. In particular, the quenched approximation, while reasonably 
successful when used in the calculation of baryon and meson spectra, appears 
to fail for the moments. The accurate calculation of the moments by lattice 
methods will presumably involve going beyond that approximation, and will 
provide a precision test of the methods used. 

We find also that the world-sheet picture of baryon structure gives useful 
insights into the moments problem, and provides a new point of view which 
could be developed further. Problems which need further study within this 
approach to QCD include the following: 

(i) establishing the connection to, and the relevance of, the chiral limit; 

(ii) the development of methods to incorporate loop effects which build in the 
extended structure of the baryons; 

(iii) and the possible usefulness of string theory methods in the calculation of 
loop effects. 

Acknowledgments 

The authors have benefited from conversations with Dr. Nora Brambilla, and ap- 
preciate her interest in this work, and her organization of the Como Conference. 
This research was supported in part by the U.S. Department of Energy under 
Grant No. DE-FG02-95ER40896, and in part by the University of Wisconsin 
Graduate School with funds from the Wisconsin Alumni Research Foundation. 

References 
References 

[1] N. Brambilla, P. Consoli, and G.M. Prosperi, Phys. Rev. D 50, 5878 (1994). 



10 



[2] J. Carlson, J. Kogut, and V.R. Pandharipande, Phys. Rev. D 27, 233 
(1983); 28, 2807 (1983). 

[3] S. Capstick and N. Isgur, Phys. Rev. D 34, 2809 (1986). 

[4] E. Jenkins, M. Luke, A.V. Manohar, and M.J. Savage, Phys. Lett. B 302, 
482 (1993). 

[5] M.A. Luty, J. March-Russell, and M. White, Phys. Rev. D 51, 2332 (1995). 
[6] J.W. Bos, D. Chang, S.C. Lee, Y.C. Lin, and H.H. Shih, |hep-ph/9602"251 . 
[7] E. Jenkins and A.V. Manohar, Phys. Lett. B 255, 558 (1991). 
[8] A. Bincer, Phys. Rev. 118, 855 (1960). 



11 



